Dynamical Critical Phenomena in three-dimensional Heisenberg Spin Glasses 



Mitsuhiro Matsumoto, Koji Hukushimajj and Hajime Takayama 
Institute for Solid State Physics, University of Tokyo, 
5-1-5 Kashiwa-no-ha, Kashiwa, Chiba 277-8581, Japan 
(Dated: February 1, 2008) 

Spin-glass (SG) and chiral-glass (CG) orderings in three dimensional (3D) Heisenberg spin glass 
with and without magnetic anisotropy are studied by using large-scale off-equilibrium Monte Carlo 
simulations. A characteristic time of relaxation, which diverges at a transition temperature in the 
thermodynamic limit, is obtained as a function of the temperature and the system size. Based on 
the finite-size scaling analysis for the relaxation time, it is found that in the isotropic Heisenberg 
spin glass, the CG phase transition occurs at a finite temperature, while the SG transition occurs 
at a lower temperature, which is compatible with zero. Our results of the anisotropic case support 
the chirality scenario for the phase transitions in the 3D Heisenberg spin glasses. 



I. INTRODUCTION 

Chirality in frustrated magnets, which provides fas- 
cinating magnetic orderings and/or novel critical phe- 
nomena, has attracted interest in recent yearsEJ. In the 
spin-glass (SG) study, the chirality-driven mechanism has 
been proposed by Kawamurao in order to explain the ori- 
gin of experimentally observed SG phase transitions. The 
chirality scenario assumes that a chiral-glass (CG) long- 
range order exists at low temperatures with preserving 
the proper rotational symmetry in a fully isotropic SG 
model. By random magnetic anisotropy even of a small 
magnitude, which always exists in real experimental situ- 
ations, a SG phase is expected to emerge as a result that 
the anisotropy mixes the two degrees of freedom, spin 
and chirality. Therefore, in the chirality-driven mecha- 
nism the SG phase transition experimentally observed in 
a class of compounds such as CuMn is essentially gov- 
erned by the CG fixed point. 

Some numerical results on equilibrium properties sup- 
port that the chirality-driven phase transitions occur in a 
three dimension^L/SI?) Edwards- Anderson (EA) Heisen- 
berg SG modem-El as well as in a 3D XY SG modeM The 
CG transition temperature and the critical exponents 
are estimated. They seem to be compatible with those 
observed experimentallyO. However, some researchers 
have argued the possibility of a finite-temperature SG 
transition^, or even simultaneous SG and CG transition 
in the 3D Heisenberg SG modeM In a 3D XY SG model, 
large scale domain-wall renormalization group calcula- 
tion at zero temperature suggests that a SG long-range 
order occurs at noruzcro temperature below the CG tran- 
sition temperaturcEJ. A similar study obtained by differ- 
ent boundary conditions also_support the existence of SG 
order at finite temperatureslij. 

The observed critical behavior associated with the CG 
transition is found to be remarkably different from that 
observed in ordinary second-order phase transitions. For 
instance, the crossing of the Binder parameter of sys- 
tems with different sizes at the transition temperature, 
which is usually observed in standard second-order tran- 
sitions, never occurs in this modcltall Furthermore, the 
previous study revealed that the CG phase of the 3D 



isotropic Heisenberg EA model had a peculiar feature 
which was consistent, with one-step replica symmetry 
breaking (RSB) phaseta. Some mean-field SG models with 
the one-step RSB exhibit dynamical singularity at a tep,- 
perature well above the static transition tcmpeeaturcE3. 
Recent simulation on a mean-field SG mode£3 shows 
that dynamical finite-size scaling holds for the dynam- 
ical transition temperature which is different from the 
static one . It is of particular interest to clarify whether 
the CG phase of the 31? Heisenberg EA model belongs 
to this class or not. 

Until now, there have been a few works on dynamical 
properties of the Heisenberg SG models 1 — .Olive et altfl 
and Yoshino and one of the authors (HT)L\3 investigated 
spin autocorrelation functions in equilibrium by Monte 
Carlo simulations. They found that, with the help of 
dynamical scaling analysis, the temperature dependence 
of the characteristic relaxation time is compatible with 
a zero-temperature SG transition. However, the corre- 
sponding CG relaxation time has not been investigated 
yet. 

In the present paper, we study dynamical properties 
associated both with spin and chiral degrees of freedom 
above the transition temperature of the 3-D Heisenberg 
EA model with and without anisotropic interactions. We 
developed an efficient method for extracting the equilib- 
rium relaxation time from non-equilibrium simulation for 
a given finite-size system and analyzed the time scale by 
using the finite-size scaling. The results of our analyses 
on the isotropic Heisenberg SG model support the first 
assumption of the chirality mechanism, namely the exis- 
tence of a finite temperature CG phase transition with- 
out a conventional SG long-range order. In particular, 
a crossover from high-tcmpcrature SG dominant regime 
to the CG dominant regime is clearly observed with de- 
creasing temperature. In addition, by introducing a finite 
random anisotropic interaction, the SG relaxation time is 
found to be significantly enhanced at low temperatures 
and in systems with large sizes, while there is no such 
drastic change in the CG dynamics. This is consistent 
with the chirality mechanism. 

Our work is distinct from the previous works in the 
following senses; (i) The CG dynamics as well as the 
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SG ones have been studied systematically both above 
and around a transition temperature, (ii) The equilib- 
rium relaxation is estimated from non-equilibrium MC 
simulations. These results may .be complementary with 
the previous equilibrium studiesuu. 

The present paper is organized as follows. In section II, 
we explain the model and the method of our Monte Carlo 
simulation. We introduce a dynamical ratio function and 
discuss its characteristic features in section [II . Finite- 



size scaling of the characteristic time is studied in section 
IV and section [v| is devoted to summary of the paper. 



II. MODEL AND SIMULATIONS 

The model studied is a classical Heisenberg model on 
a simple cubic lattice, defined by the Hamiltonian, 

H(S) = -J2JiA-Sj- £ DgSfS], (1) 

where Si = (Sf, S? , Sf ) is a three-component unit vector, 
and the sum runs over all nearest-neighbor pairs with the 
total number of sites TV = LxLxL. The nearest neighbor 
couplings Jij are random bimodal variables which take 
± J with equal probability. The anisotropic interactions 
are also independent random variables obeying a 
box distribution with its range [-D : D]. We assume 
that the anisotropic term is symmetric, D?f = D"? . A 
quite similar model has been studied by MC simulationEl 
with some values of strength D, where the SG phase 
at low temperatures has been found around and above 
D/J^ 0.05. We study the model system both with 
the isotropic case (D/J = 0) and an anisotropic case 
(D/ J = 0.05). 

We employ a standard heat-bath Monte Carlo method 
with two sub-lattice flipaU. One Monte Carlo step 
(MCS) is defined as TV spin trials. We perform MC pro- 
cedure as follows; firstly spins are set as a random config- 
uration, namely they are instantaneously quenched from 
the high temperature limit. The spin configuration is up- 
dated at a working temperature T during a waiting time 
t w , which is from 10° to 10 5 . Subsequently, additional 
MC steps are performed for measurements. The lattice 
sizes L studied are 4, 8, 16 and 32. In our simulations, we 
perform eight independent runs for each identical sample 
to take an average over initial conditions and different 
thermal noises. A typical number of samples averaged 
over is 40-960 depending on the system size. 

In order to study non-equilibrium dynamics of the 
model, we calculate the spin autocorrelation function de- 
fined by 



Cq(t Wl t) — — ^ ^ [(Sj(t w ) ■ Si(t y 



t))h 



(2) 



where (• • ■) represents an average over initial conditions 
and [• • •] over the bond disorder. In the long time limit 
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FIG. 1: Time dependence of the spin-g[ass (upper) and chiral- 
g[ass (lower) autocorre[ation functions with different waiting 
times in the 3D isotropic Heisenberg SG modei at T/J — 0.20. 
The system size is fixed, L = 16. 



(t — » oo) after taking the equilibrium limit (t w — > oo) the 
autocorrelation, Eq. (||), converges to the SG order pa- 
rameter originally proposed by Edwards and AndcrsonliZI. 
We also calculate the corresponding chirality autocorre- 
lation function, which allows us to detect a possible CG 
transition from the dynamical point of view. It is also 
defined as an overlap of the local chirality 



Xi(tw+t))], 



(3) 



where Xiu is the local chirality at the ith site and in the 
/ith direction defined by a scalar as 



Xi 



S 



[Si X Si — e fl ), 



(4) 



with e M (/i = x,y,z) being a unit lattice vector along 
the /j, axis. The CG autocorrelation function in off- 
equilibrium has not been investigated extensively yet, ex- 
cept for Ref. |[ 

In Fig. |l|, we show the SG and CG autocorrelation 
functions with the different waiting times t w against time 
t. In non-equilibrium dynamics, temporal correlation 
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FIG. 2: Time dependence of the spin glass (upper) and chiral- 
glass (lower) autocorrelation functions with different sizes and 
the fixed waiting time t w = 10 s in the 3D isotropic Heisenberg 
SG model at T/J = 0.20. 



functions explicitly depend on not only t but also t w , 
which we call the breaking of time translational invari- 
ance. The longer is the waiting time, the slower be- 
comes the relaxation. The curves with shorter waiting 
times start to deviate from an equilibrium function at the 
shorter time scales. Such a time scale, below which the 
time translational invariance holds effectively, is roughly 
proportional to the waiting time t w . 

In Fig. |^ we present the correlation functions with dif- 
ferent system sizes for a fixed waiting time t w = 10 5 . One 
can clearly see that the relaxation curves start to deviate 
from that of the thermodynamic limit. This implies that 
finite-size effects yield another characteristic relaxation 
time. The latter is related to the critical slowing down 
which we study in the present work. 

Generally speaking, equilibrium MC simulations for 
SG models suffer from extremely slow dynamics and equi- 
libration is hardly realized in large systems at low tem- 
peratures. One of the advantages of the non-equilibrium 
simulation is to practically avoid such a very slow equi- 
libration process. In the present work, we focus our at- 
tention to the longest relaxation time of finite-size sys- 



tems which is expected to diverge with the system size at 
and below the transition temperature. Our strategy of 
the non-equilibrium simulation is to take the equilibrium 
limit, i.e., t w — > oo, after extracting the characteristic re- 
laxation time under off-equilibrium in finite systems and 
then take the thermodynamic limit L — > oo. 



III. DYNAMICAL RATIO FUNCTION AND 
RELAXATION TIME 

According to the dynamic scaling ansatzEH, the auto- 
correlation function in equilibrium near a second-order 
phase transition has the form 



C(t w = oo,t w + t) ~r x f(t/T) 



(5) 



where r is the relaxation time at a given temperature and 
A a critical exponent. This form represents a crossover 
around r from the short-time critical behavior, to 
the long-time relaxation described by f(t/r). When the 
temperature approaches to the critical point T c , the re- 
laxation time t diverges as r(T) ~ \T — T c \~ zu . Thus, 
the autocorrelation function follows a power law in t at 
T c . 

One may consider that the critical temperature as well 
as the critical exponents associated with the transition 
can be easily obtained by the relaxation time near the 
critical temperature. There are, however, several diffi- 
culties in obtaining the relaxation time by MC simula- 
tions. One difficulty is due to extremely slow dynamics 
of SG systems as mentioned above. Another difficulty is 
that the method of obtaining the relaxation time requires 
accurate estimation of a long-time tail of the autocorre- 
lation function. One way is to integrate C{t) over time. 
A contribution to the integrated relaxation time is dom- 
inated by the long-time tail where C (t) simulated is very 
small and largely fluctuating. Another method would be 
a fitting of C(t) to a relaxation function with the scaled 
argument t/r. It is, however, difficult to know the func- 
tional form explicitly, which is often non-exponential as 
seen in Fig. [l]. 

In order to avoid the above mentioned difficulties, 
Bhatt and Youngll3 (hereafter referred to as BY) have 
introduced a dimensionless dynamic correlation function 



W X/i Si(t w ) ■ Si(t w + t) 



_ . (6) 
The prefactor in Eq. (|£j) involving a power of t is canceled 
out by taking the ratio of moments. A similar dimension- 
less dynamical, function has also been studied in simple 
ferromagnetscJ'El. The corresponding CG ratio function 
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FIG. 3: Double logarithmic plot of the spin-glass (upper) 
and chiral-glass (lower) autocorrelation functions and ratio 
functions of the isotropic 3D Heisenberg spin glass model at 
T/J = 0.30. 
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FIG. 4: Time dependence of the SG (upper) and CG (lower) 
ratio functions with different waiting times in the 3D isotropic 
Heisenberg SG model at T/J = 0.20. 



is defined as 



(7) 

Because these ratio functions are dimcnsionless such as 
the Binder parameter in the static case, the dynamical 
scaling form is expected to be given by a function of t/r 

as 



R(t) ~ R(t/r), 



(8) 



where R(x) is a universal scaling function. BYE2I fixed 
the waiting time t w in Eq. (^) to t because it was be- 
lieved that equilibration of the autocorrelation function 
on a given time scale t needed the waiting time t w of 
the same time scale as t. They studied the dimensionless 
function of the short range Ising SG and the mean-field 
Sherrington-Kirkpatrick models just at the known criti- 
cal temperature and successfully obtained the dynamical 
critical exponent z using the finite-size scaling for the 
relaxation times t(L, T c ) ~ L z . 



In the present study, we systematically investigate the 
ratio function both for the SG and CG autocorrelation 
functions as bivariate functions of t and t w . Instead of 
taking the equilibrium limit of the ratio function or fix- 
ing the time scale t w ~ t studied previously, we try to 
take the equilibrium limit (t w — > oo) of the relaxation 
time r(t,t w ) obtained from the ratio function for several 
waiting times. 

Equations (^J) and (Q) are the ratio of an odd moment 
of the autocorrelation function to an even one. Conse- 
quently, the ratio functions for a given t w decay from 1 
at t = to zero as t goes to infinity and they are sen- 
sitive to flips or rotations of the entire system, which 
are the longest mode of relaxation. Therefore we expect 
that the ratio functions pick up only slowest relaxation 
modes from the whole modes of C(t). A typical exam- 
ple of the ratio function, as well as C(t), is presented for 
the SG and CG autocorrelations in Fig. ||. As expected, 
the ratio functions are unity in the relatively short-time 
regime, implying that the denominator and the numera- 
tor coincide with each other. For the longer-time regime, 
as shown in the figure, the ratio functions of spin and chi- 
rality sector coincide with the tail of the corresponding 
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FIG. 5: Scaling plot of the SG(upper) and CG(lower) ratio 
functions of the 3D Heisenberg SG model with various tem- 
peratures, waiting times and system sizes. 

autocorrelation function multiplied by a certain constant, 
which represents a statistical weight of the slowest modes 
in C{t). 

When the ratio function R(t) extracts only one relax- 
ation mode corresponding to the longest mode in the 
system, it is expected to show a universal scaling form 
and to be scaled as 

R(t w ,t w +t;T,L) ~ R(t/r(t w ;T,L)). (9) 

Here we assume that the ratio functions depend on t w , T 
and L only through the relaxation time. This may be 
an extension of the dynamical scaling (||) to an off- 
equilibrium situation, which is to be checked. We try to 
scale the observed ratio functions into a universal scaling 
function by appropriately choosingthe time constant r 
depending on t w , T and L. Figure g presents a result of 
the scaled ratio functions for D/J = 0. In the figure, 
one sees that both the SG and CG ratio functions with 
various temperatures, system sizes and waiting times are 
almost scaled as a common universal function. Note that 
this finding allows us to compare the absolute values of 
the SG and CG relaxation time. Using the scaling of the 
ratio functions in this way, we obtain the characteristic 



FIG. 6: Temperature dependence of the SG (upper) and CG 
(lower) relaxation times in equilibrium for the isotropic case 
(open marks) and anisotropic case with D/J — 0.05 (filled 
marks) . 

relaxation times r(t w ;T,L) for each set of parameters 
t w ,T and L. We then let t w to go to infinity to obtain 
the characteristic relaxation times in equilibrium. We as- 
sume that the t^-dependence of the relaxation times is 
described as 

T(t w ;T,L)=T(t w = oo;T,L)-at-P, (10) 

with a,/? being some constants and t(T,L) = r(t w = 
oo;T, L) the equilibrium relaxation time. 

Temperature dependences of t(T, L) thus extracted are 
shown in Fig. |6| for the spin (upper) and chiral (lower) 
degrees of freedom. We first note that the CG relaxation 
time is almost independent of the strength of the mag- 
netic anisotropy D. This is natural since the anisotropy 
breaks the rotational symmetry, but not the reflection 
symmetry which is spontaneously broken by the CG 
phase transition. On the other hand, as seen in Fig. [(], 
the random anisotropy changes drastically the behavior 
of the SG relaxation time. The SG relaxation time of 
D/J = 0.05 increases rapidly around T/J = 0.25 with 
decreasing the temperature. Such a drastic change of the 
behavior implies that the anisotropy mixes the chirality 
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with spin degrees of freedom as predicted in the chirality 
scenario. This is the first observation of the effect of the 
anisotropy for the SG relaxation times. 

We can show another evidence that the chirality domi- 
nates the phase transition in the isotropic 3D Heisenberg 
SG model. In Fig. [?], the time evolution of SG (bro- 
ken line) and CG (solid line) ratio functions at T/J = 
0.40,0.30 and 0.20 with L = 16 and t w = 10 5 are shown 
simultaneously. The characteristic relaxation time of the 
CG ratio function becomes larger with decreasing tem- 
perature much more rapidly than that of the SG ratio 
function does and it finally exceeds the SG relaxation 
time at T/J = 0.20. This fact strongly suggests that the 
CG degrees of freedom dominate large-scale and long- 
time phenomena at low temperatures, while the SG one 
dominates at relatively high temperatures. This type 
of crossover is observed at the first time in the present 
work. We can also see such a crossover in the resulting 



equilibrium relaxation time. Figure [3| shows the same 
data as shown in Fig. ^ but only for the isotropic case 
D/J = 0. A characteristic temperature at which the 
two relaxation times cross each other shifts towards the 
high-temperature side as the system size increases. We 
emphasize that this crossover is not due to an apparent 
finite-size effect. 

IV. FINITE-SIZE SCALING 

In the present section we discuss the finite-size-scaling 
analysis of the equilibrium relaxation time obtained in 
the previous section. According to a standard dynamical 
finite-size-scaling analysis, the scaling form of t(T, L) is 
given by 

r(T, L) ~ L z t ((T — X' C )L 1/V ^ , (11) 

where z is the dynamical critical exponent and v the 
correlation-length exponent. The scaling function be- 
comes a constant just at the critical temperature, lead- 
ing the .finite-size scaling form, t(T c , L) ~ L z , previously 
studicdEl We discuss asymptotic behavior of the scaling 
function t(x). In a large-system-size limit above T c , i.e., 
i>1, the scaling of Eq. (jll]) should recover an ordinary 
scaling form, 

t(T,L = oq)~(T-T c )-"'. (12) 

This leads that the scaling function behaves as 

t{x) ~ x~ zu (13) 

at x 3> 1. Assuming that the transition takes place at 
a finite temperature, the opposite limit, i.e. —x 3> 1, 
can be also argued in the isotropic case where the sys- 
tem would entirely behave like a rigid magnet. For the 
SG ordering in D/J = 0, such a system rotates glob- 
ally as coherent diffusive motion with the relaxation time 
td proportional to the bulk volume. This implies that 
the asymptotic scaling form of the SG relaxation time in 
D/J = is 

t{x) ~ \x\^ d ~ z)v , (14) 

in the limit —x 3> 1 with d being the dimensionality. The 
asymptotic scaling function for the CG relaxation time, 
however, would follow an exponential activated type be- 
cause of the discrete nature of chirality. 

The situation is complicated in a system with T c = 0. 
In one case where the conventional form ( pT[ ) holds with 
T c = as in the standard second-order phase tran- 
sition, the scaling function in the asymptotic limit at 
x <C 1, t(x = 0), becomes .-a, constant. On the other 
hand, Yoshino and TakayamaliJ proposed a modified scal- 
ing form which was a combination of the standard form 
( pd| ) with T c = and the diffusion relaxation time, 
t d ~ L d /T. This form is compatible with (||) and (0) 
in the limit x •< 1 if (d — z)v = —1 holds. In this case, 
the scaling exponents, z and v, arc no longer independent 
parameters. 
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FIG. 9: Finite-size scaling plot of the SG relaxation time in 
the isotropic 3D Heisenberg SG model. The parameters of 
the scaling are as follows : Tsg = 0.0(1), v = 1.4(2) and 
2 = 4.4(4). The solid line represents the expected asymptotic 
behavior for large scaling regime. 



A. Isotropic case: D/J = 

The result of the finite-size scaling analysis mentioned 
above is demonstrated for the SG relaxation times of the 
isotropic Heisenberg SG model in Fig. |[ As can be seen 
from the figure, the data points are found to collapse into 
the universal scaling function. From the scaling we esti- 
mate T SG = 0.0(1), v = 1.4(2) and z = 4.4(4). Our result 
for the SG transition temperature is roughly .CQaaigteai: 
with the zero-temperature phase transitiorolijIijEaEi 
although we could not completely exclude the possibility 
of the SG phase transition at a very low but finite tem- 
perature. Assuming the zero-temperature transition, the 
scaling results suggest that relaxation times diverge with 
a power law of T as T goes to 0, but not an exponential 
divergence which occurs at the lower critical dimensions 
such as the two-dimensional Ising SG model. 

We expect that the scaling function converges to a fi- 
nite value in the limit x — > as in the standard second- 
order phase transition, though we don't have enough data 
to confirm this convergence. The scaling plot may be 
consistent with a modified scaling from x > 1 to x < 1 
proposed by Yoshino and Takayamatj, but then temper- 
ature dependence of the relaxation time at low tempera- 
tures, t ~ j'(<i-2)i> w jth z and v obtained above does not 
coincide with that expected from a simple diffusive relax- 
ation time proportional to 1/T. A further investigation 
is needed to clarify this point. 

The present estimate for v is compatible with the. pre- 
vious ones in other works; v = 1.54(19)t3 and 2.0(2)0 by a 
domain- wall renormalization group calculation at T = 0, 
v = 1.35(5) by a^HC simulatior£3, v ~ 1.6 by an eqp,- 
librium dynamicsEj and v ~ 1.14 by a MC simulationO. 

We have also analyzed the CG relaxation time by a 
similar scaling. The best scaling shown in Fig. HO is 



FIG. 10: Finite-size scaling plot of the CG relaxation time 
in the isotropic 3D Heisenberg SG model. The parameters 
of the scaling are as follows; Tcg = 0.22(2), v = 1.1(3) and 
z = 4.1(4). 




1 10 100 1000 10000 100000 

f[MCS] 

FIG. 11: Time dependence of C x /R x for the CG at T/J = 
0.20. 



obtained by the following estimates: Tcg = 0.22(2), 
v = 1.1(3) and z = 4.1(4). The estimate for Tqg ^r^fFr 
sistent with the recent equilibrium MC calculationsoEa. 
The present estimate for v is rather close to the previous 
MC estimate for the Gaussian distribution, v ~ 1.2a 

At the critical temperature, the autocorrelation func- 
tion C x and also C' x /R x follow the algebraic decay (|5|) 
with the exponent A = (3/zv. In Fig. [Tljwe show C x /R x 
as a function of time at T/J = 0.20 close to the CG 
transition temperature. The data still depend on both 
the system size and the waiting time. As a preliminary 
study, we estimate the exponent A = 0.23(5) from the 
slope of the plot, where the numerical error comes from 
the uncertainty of the critical temperature. This value 
combined with z and v by the finite-size scaling leads to 
another exponent, f3 ~ 1.0, which is consistent with the 
previous workcl. 

As mentioned in the introduction, the dynamical phase 
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transition which is separated from the static one can ex- 
hibit at a finite temperature in some one-step RSB sys- 
tems. It depends on whether the order parameter contin- 
uously appears or discontinuously jumps up at the static 
transition as the temperature decreases. Examples of the 
discontinuous one-step RSB transition are the mean-field 
p— spin glass with p > 2 and the mean-field g-state Potts 
glass with q > 4, while those of the continuous one-step 
RSB transition are the mean-field g-state Potts glass with 
2.8 < q < 4. Our estimation of the CG transition tem- 
perature by the dynamical finite-size scaling is consistent 
with the static transition temperature. This fact means 
that the CG phase belongs to the latter class of the one- 
step RSB. 

B. Anisotropic case : D/J = 0.05 

The finite-size scaling plots of the CG relaxation times 
in the anisotropic 3D Heisenbcrg SG model is shown in 
Fig. |l2|. In the analysis we regard the CG critical tem- 
perature as a unique adjustable parameter, and obtain 
Tcg = 0.24(2). As seen in Fig. [j^, the scaling analysis 
works well for both D/J = 0.0 and 0.05 with the com- 
mon z and v. This implies that the critical exponents 
associated with the CG critical phenomena as well as the 
scaling function arc independent of the anisotropic term. 
This finding of our analysis is consistent with the expec- 
tation that the anisotropic interactions are expected not 
to affect the CG symmetry. 

We believe that the SG phase transition also occurs 
in the anisotropic case at the same temperature as the 
CG transition temperature. In fact, we have found in 
Fig. |^ that the anisotropic term significantly enhances the 
SG relaxation time at low temperatures. Unfortunately, 
however, such enhanced relaxation time appears only in 
larger sizes and at lower temperatures. Therefore, these 
data points are not enough to test the finite-size scaling. 



V. DISCUSSIONS 

We have studied dynamical critical phenomena in the 
3D Heisenberg SG models using Monte Carlo simula- 
tions. Our results support the chirality scenario for re- 
alistic SG transitions. It is confirmed that in the fully 
isotropic model at low temperatures, the CG relaxation 
time exceeds the SG relaxation time. This result strongly 
suggests that long-time or large scale behavior is domi- 
nated by the CG ordering. The finite-size-scaling anal- 
ysis of the relaxation time supports that the CG phase 
transition with the diverging relaxation time occurs at a 
finite temperature, while the SG phase transition takes 
place at a lower temperature, consiai^pt-atiih even at zero 
temperature as previously studicdEjOlijEjO. In other 
words, the spontaneous breaking of discrete Z2 symme- 
try likely occurs with preserving the proper rotational 
symmetry SO(3). The SG and CG transition tempera- 




0.1 1 10 



\T-T CG \L UV 

FIG. 12: Finite-size scaling plot of the CG relaxation time 
in the 3D Heisenberg SG model with D/J = 0.05. The 
parameter of the scaling is the CG transition temperature 
T CG = 0.24(2). The scaling data of D/J = shown in Fig. [k] 
are also plotted with the cross mark. 



tures of our result are significantly separated with each 
other within numerical accuracy. This is the first support 
for the basic assumption of the chirality scenario from a 
dynamical point of view. 

We have also studied an anisotropic 3D Heisenberg SG 
model in the similar way. It is found that by introducing 
the weak random anisotropy, the SG relaxation times 
in the large system sizes increase at low temperatures, 
similar to that observed in the CG relaxation time, al- 
though a finite-temperature SG phase transition cannot 
be directly determined by the finite-size-scaling analysis. 
Our finite-size-scaling analysis of the CG relaxation time 
strongly suggests, on the other hand, that the CG phase 
transition with the anisotropic interactions belongs to the 
same universality class as the isotropic case. When the 
SG phase transition in the anisotropic case is intrinsically 
dominated by the SG order parameter but not the CG or- 
der parameter, the observed CG critical behavior should 
be affected by the SG transition and it should be differ- 
ent from the fully isotropic case. Therefore, from this 
point of view, our result is consistent with the chirality 
mechanism where the phase transition is dominated by 
the CG order parameter. In order to further clarify the 
nature of phase transition in the anisotropic case, how- 
ever, directly measurements of the SG phase transition 
and its critical exponents are required. 

A standard scenariccj based on the zero-temperature 
SG transition in the isotropic point says that the SG 
transition continuously approaches to zero in the zero 
anisotropy limit as Tsg(-D) oc D 1 / 4 ; see Fig. |l3|(a). Be- 
cause the anisotropy reduces the symmetry of Hamilto- 
nian, the universality class changes from the isotropic SG 
to the anisotropic SG. The latter is expected to be in the 
class of the Ising SG n On the other hand, according to the 
chirality mechanismel, the random magnetic anisotropy D 
plays a role of mixing the CG long-range order with the 
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D/J D/J 

FIG. 13: A schematic picture of a phase diagram in the stan- 
dard scenario (left) and the chirality scenario (right). 



limit and the universality class along the phase bound- 
ary. Our findings that the CG universality class is not 
affected by the anisotropy D and that Tsg < Tcg m the 
isotropic system prefer the chirality scenario. 

In summary, spin-glass and chiral-glass orderings in 
three-dimensional Heisenberg spin glass model are stud- 
ied from the point of view of dynamics by off-equilibrium 
Monte Carlo simulations. Our results support the chiral- 
ity scenario for the SG phase transition in three dimen- 
sions. 



SG one. The SG transition temperature in the ideal zero 
anisotropy limit is predicted to coincide with Tcg {D = 0) 
but not Tsg(D = 0), as shown in Fig. |l^(b). This may 
be true wherever Tcg > Tsg at the isotropic limit even 
with Tsg(D = 0) > 0. Once the CG has a long-range 
order, even the small anisotropy D coherently breaks the 
rotational symmetry of the system. Consequently, the 
SG critical behavior observed in the anisotropic model is 
governed by the fixed point of the CG transition of the 
isotropic model. The two scenarios can in principle be 
distinguished by the phase boundary near the isotropic 
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